Level statistics and localization in a 2D quantum percolation problem 
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A two dimensional model for quantum percolation with variable tunneling range is studied. For 
this purpose the Lifshitz model is considered where the disorder enters the Hamiltonian via the 
nondiagonal elements. We employ a numerical method to analyze the level statistics of this model. It 
turns out that the level repulsion is strongest around the percolation threshold. As we go away from 
the maximum level repulsion a crossover from a GOE type behavior to a Poisson like distribution is 
indicated. The localization properties are calculated by using the sensitivity to boundary conditions 
and we find a strong crossover from localized to delocalized states. 
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The statistics of energy levels in complicated quan- 
tum systems has been a subject of research for several 
decades. It started with the study of energy levels in 
atomic nuclei 0| , then the statistics of electronic states in 
atoms was investigated |Q , and more recently the statis- 
tics of electrons in quantum dots (also known as "ar- 
tificial" atoms) H] . A surprising result of most of these 
approaches was that the statistics of energy levels is quite 
universal regardless of the specific system: systems can 
be classified according to their symmetry properties as 
orthogonal, unitary and symplectic. These universality 
classes represent strong correlations between the energy 
levels due to level repulsion. This is indicated, e.g., by 
the level spacing distribution P{s) which goes like 
for small s. The exponent (3 (=1,2,4) characterizes the 
universality class Q]. In contrast to this repulsive be- 
havior (Wigner-Dyson distribution) the assumption of 
statistically independent energy levels would lead to a 
Poisson distribution P{s) = e~*. The correlation in nu- 
clei or atoms is so strong because their corresponding 
states have usually a large overlap, except perhaps for 
the ground state. The situation is less clear if we consider 
a macroscopic system of atoms like in solid state physics 
where disorder can be present. Depending on the latter 
there are extended electronic states and also states which 
are localized in space due to Anderson localization |^ . In 
terms of level statistics the localized states are expected 
to obey a Poisson distribution, whereas the extended 
states are expected to have strong level repulsion char- 
acterized by a Wigner-Dyson distribution. Since there 
is a metal-insulator transition from extended to local- 
ized states driven by disorder (Anderson transition) , it is 
natural to study this in terms of the statistics of energy 
levels. 

Approaches to the role of level statistics for the An- 
derson transition |^ were presented by several authors, 
and a model for quantum percolation for nearest neigh- 
bor transfer was investigated in Ref . (D . Recently it was 
found that the divergence of the localization length 
^ at the metal-insulator transition leads to a deviation 
from the Wigner-Dyson statistics for £^ > L, L being the 



system size. In particular it was found that the decay 
of the level spacing density is weaker than that of the 
Wigner-Dyson statistics. This is also found in our inves- 
tigation. 

The purpose of this paper is to study a quantum per- 
colation model where the transfer is not only between 
nearest neighbors but where the transfer rate decays ex- 
ponentially with distance. Moreover, in contrast to An- 
derson's model for localization we study disorder in the 
off-diagonal part of the Hamiltonian. This is also known 
as Lifshitz type of disorder. This model is motivated by 
physical systems. One example is a two-dimensional ar- 
ray of quantum dots |^ . Our model is also motivated by 
the analogous picture of variable range hopping in solids 
A third example for this model is the class of low 
doped high-Tc cuprates. Here the charge carriers are 
holes in two-dimensional Cu02 layers. It seems that dis- 
order and phase separation play an essential role in these 
systems [|lO|,^ . The two-dimensional copper oxide plane 
is separated into hole rich conducting and magnetically 
correlated (insulating) areas. A possible origin for phase 
separation are polaronic states that are discussed in Refs. 
|l2| , ^ . Therefore our model can describe a transition or 
at least a crossover from strongly insulating (localized) 
states to states with infinite or at least very large local- 
ization length. This picture can be applied to the physics 
of the normal state in the low-doped high-Tc materials. 
Indeed, conductivity measurements by Chen et al. p^ ] 
of low doped La2-2;Sr2;Cu04 with x « 0.002 have shown 
that for temperatures below 50 K the transport proper- 
ties are governed by a hopping type conduction. Also 
earlier measurements by Keimer et al. for a sample 
with doping concentrations of x w 0.04 found a conduc- 
tivity of hopping type near localization below 20 K. 

Our model for quantum percolation corresponds to the 
Hamiltonian 



(1) 



with the following off-diagonal (hopping) matrix ele- 
ments 
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te^p{-a{rij - tq)) for r,y > tq 
t for rij < ro 
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The lattice sites i and j can be randomly occupied with 
quasiparticles in Wannier states, leading to random hop- 
ping elements tij . This type of randomness is also known 
as Lifshitz type disorder The exponential decay of 
the localized wave functions leads to an exponentially 
decay of the hopping rate with distance on the inverse 
decay length a. The spatial extension of the localized 
states, e.g., given by the size of a polaron, is expressed 
with tq. a hopping matrix element tij is non-zero only 
if the sites i and j of the 2D lattice are both occupied by 
localized states. 

An advantage of the long range hopping of our model, 
at least for small enough a, is the fact that the density 
of states is smoothed out in contrast to the sharp peaks 
found for nearest neighbor transfer Q . The smooth den- 
sity of states is easier to analyze with the methods of 
random matrix theory. The density of states is shown in 
fig. H It shows a broad peak near the lowest eigenvalue. 
This peak is due to the two-dimensional nature of the 
system. For infinite a and all lattice sites occupied only 
a nearest neighbor transfer remains. In this case the den- 
sity of states is the elliptical integral with the logarithmic 
singularity at the center. For finite a also next nearest 
neighbor and further transfers are included. This shifts 
the peak in the density of states to the lower band edge. 
The density of states drawn in fig. ^ shows a remnant of 
this peak. 

The numerical calculation is performed as follows. The 
A'' (typically N = 400, N < L^) localized states are 
randomly chosen with probability c (c = N/L'^) on an 
L X L square lattice with lattice constant a. In this pro- 
cedure periodic boundary conditions are used. The coor- 
dinates of the localized states are distributed randomly 
while multiple occupation is prohibited. 

For example, for ro = a the classical (bond) percola- 
tion threshold is near c « 0.5. For any combination of 
pairs the off-diagonal elements have to be computed for 
the Hamiltonian (|l|) with the hopping element (^). The 
resulting matrix is diagonalizcd numerically using stan- 
dard orthogonal decomposition methods. It is important 
to notice that in contrast to the corresponding matrix of 
the Anderson model we do not obtain a sparsely occu- 
pied matrices. This requires more numerical effort and 
leads to a limitation of the matrix size we can diagonalize 
to matrices not larger than 400 x 400. As a result the 
distribution curves fluctuate stronger than in the case of 
nearest neighbor hopping models, where the matrices can 
be significantly larger. 

The level spacing distribution P{s) of our model is an- 
alyzed and compared with the Poisson distribution and 
with the distribution of the Gaussian orthogonal ensem- 
ble (GOE, P^{s) = Trfe"'^^). The choice of the GOE is 
due to the fact that our Hamiltonian obeys time reversal 
symmetry. As one can see in fig. ^ the level-statistics 
for the quantum percolation regime does not follow the 



GOE regime. In particular, for s > 2 the distribution 
decays slower than P'^(s). This is in agreement with the 
prediction by Aronov et al. ^ for the situation near the 
Anderson type metal insulator transition. It can be in- 
terpreted as the domination of the statistics for larger 
level spacings by weakly overlapping states. According 
to our results this is not only a feature near the metal 
insulator transition but it is present in the whole dop- 
ing range. In general the level repulsion of the quantum 
percolation model is weaker than the one for the GOE. 
In order to investigate this behavior in more detail we 
analyze the A3 statistics of the eigenvalue spectra. The 
latter is defined for a different number n of levels as Q 

A3(n) = - f\st{x') ~ Ax' - Bfdx' . (3) 
nA,BjQ 

Here St(x') is the staircase function. The result is shown 
in fig. 0. The A3 statistics indicates the following be- 
havior for the doping dependence of the system: for low 
doping, far below the classical percolation threshold, the 
system shows only weak level repulsion and a tendency 
towards Poisson statistics. This is expected from com- 
mon arguments because well-separated localized states 
are almost independently distributed, leading to a Pois- 
son distribution. For moderate doping, in the vicinity 
of the classical percolation threshold, the level repulsion 
increases and the system shows a tendency towards the 
Wigner statistics. This behavior indicates the beginning 
of the formation of overlaps between the states. As a re- 
sult the eigenvalues experience level repulsion. However, 
this tendency is reversed, when the doping concentration 
is further increased above the percolation threshold. A 
possible explanation is a tendency towards uncorrelated 
fc-space states of the fully doped (pure two-dimensional) 
system. 

To investigate the transition between spatially local- 
ized states and extended states the sensitivity of the 
eigenvalues with respect to a change of the boundary 
conditions is considered jl7|. The Peierls substitution 
of the hopping matrix elements t — > ie"^^^ is used to 
vary continuously the boundary conditions in the Hamil- 
tonian. Expanding the exponential function enables us 
to use perturbation theory |18|| 



H 



H + H, 



= ^^{itijixi - x.j)(j)c\cj + h.c), 



(4) 



(5) 



where Xj is the x-coordinate of site j . Since is purely 
imaginary the eigenvalues of the Hermitean Hamiltonian 
are affected only in second order perturbation theory 
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For the numerical calculation the average is taken over an 
ensemble of (typically 30) matrices. N runs from M — 10 
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to M + 10 since the nearest energy levels contribute 
mainly to Ai?. A_Em can be identified with the conduc- 
tivity via the Kubo-Greenwood formula |l9| . Eqn. is 
also known as the Thouless formula for the conductivity. 

A numerical investigation can only give information 
about localization lengths smaller or comparable with the 
system size . The transition from localized to delocal- 
ized states, shown in fig. ^, may indicate a real transition 
to delocalized states in the infinite system. At least it will 
indicate a transition from exponentially to algebraically 
decaying states. Furthermore, fig. ^ shows that with in- 
creasing a the crossover from localized to extended states 
gets more abrupt while the absolute value of the energy 
shift decreases, a^^ , the characteristic length scale of 
the hopping processes, is always much smaller than the 
system size. 

In conclusion, we find a clear indication of a qualita- 
tive change of the system in the level statistics as we go 
through the percolation threshold as shown in figs. 
This effect depends on the strength of doping c as well 
as on the range of the transfer a. To explain the on- 
set of delocalization for the normal state of the high-T^ 
cuprates with our model the diameter of the polaronic 
states should be chosen as rg ~ 4a. This size is mo- 
tivated by spatial inhomogeneities seen in experimental 
observations (see e.g. inelastic neutron scattering data of 
ErBa2Cu30j; by Mesot et al. which are interpreted 
with similar cluster sizes). Such polaronic states will give 
a percolation threshold for a doping concentration of c w 
0.05. Therefore, the transition from localized to delocal- 
ized states will occur near this concentration if ~ 4a. 
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Figure captions 

FIG. 1. For different concentrations, below and above the 
classical percolation threshold, the level spacing distribution 
is drawn. It is compared with the Wigner and the Poisson 
distribution. Note that for s > 2 the distribution decays 
slower than predicted by the Wigner distribution. The pa- 
rameters for this plot are a — 0.5 /a with the lattice constant 
a and ro = a For comparison the inset shows the result for 
the GOE ensemble of the same size. 
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FIG. 2. For the same concentrations and parameters as in 
fig. 1 the results of the calculated A3 statistic is plotted. The 
transition from a more Poisson like to a more Wigner like and 
back to a Poisson like distribution can be clearly seen in this 
plot. The level repulsion is strongest near the classical perco- 
lation threshold. Again for comparison the GOE ensemble is 
shown in the inset. 

FIG. 3. The density of states for two different concentra- 
tions is drawn. The increase of the density of states towards 
small energy values is due to the two-dimensionality of the 
system. Since a = 0.5/a as in fig. 1 and in fig. 2 was chosen 
no peaks from parts separated from the percolation backbone 
occur. 
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FIG. 4. For different values of a (in units of the inverse 
lattice constant) the sensitivity to the boundary conditions is 

shown in dependence of the doping concentration c. A tran- 
sition from localized to delocalized states can be seen which 
for larger a gets more abrupt as can be seen in the inset. For 
this calculation ro is again ro = a 
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